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We prove multiplicativity of maximal output p norm of classical noise channels and thermal noise channels 
of arbitrary modes for all p > 1 under the assumption that the input signal states are Gaussian states. As a direct 
consequence, we also show the additivity of the minimal output entropy and that of the energy-constrained 
Holevo capacity for those Gaussian channels under Gaussian inputs. To the best of our knowledge, newly 
discovered majorization relation on symplectic eigenvalues, which is also of independent interest, plays a central 
role in the proof. 
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I. INTRODUCTION 

One of the goals of quantum information theory is to clarify 
the ultimate capability of information processing harnessed 
by using quantum mechanics [1, 2J. The celebrated Holevo- 
Schumacher- Westmoreland theorem 01 ^ gives us a formal 
basis to determine the ultimate transmission rate of classical 
information encoded in quantum states transmitted through a 
quantum channel. Yet, an important question is still unan- 
swered in terms of the classical capacity of quantum chan- 
nels. It is the additivity question; Do the entangled inputs 
over several invocations of quantum channels improve the 
classical capacity of quantum channels? Despite many ef- 
forts devoted to the additivity problems of quantum channels, 
the additivity properties have been proven for a few examples, 
such as entanglement breaking channels ^5J, unital qubit chan- 
nels 01, depolarizing channels 01, and contravariant channels 
iP] . Surprisingly, the additivity problems of quantum channels 
have been shown to be equivalent to the seemingly unrelated 
additivity problems of quantum entanglement, i.e., the addi- 
tivity and the strong superadditivity of entanglement of for- 
mation l9, fTolfTllfl2il . All of them are not completely solved 
and are now major concerns in quantum information and the 
quantum entanglement theories. 

As for continuous-variable quantum systems, in spite of in- 
tensive research |13, 14, 15], only lossy channels have been 
proven to be additive Id 61 . The additivity problems may be 
much more intractable for continuous-variable quantum sys- 
tems. The natural question is, therefore, what can we say 
about the additivity properties of Gaussian quantum channels 
if we restrict the input signal states to be Gaussian states? This 
question has its own significance. One rationale is that the 
Gaussian channels correspond to the so-called Gaussian oper- 
ations that can be implemented by current experimental tech- 
niques, such as bemsplitters, phase shifters, squeezers, and 
homodyne measurements. Another is the mathematical sim- 
plicity; Gaussian operations on Gaussian states are completely 
characterized by finite dimensional matrices and vectors, al- 
though the underlying Hilbert space is infinite dimensional. 
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Due to their mathematical simplicity, the additivity problems 
of Gaussian channels under Gaussian inputs provide a poten- 
tial firm step towards answering the additivity questions. 

Serafini et al. fv^ formulated the multiplicativity problems 
of the purity at the output of Gaussian channels measured by 
the Schatten p norm under the assumption that the input sig- 
nal states were Gaussian states. In this paper, we extend their 
formalism to the additivity problems of minimal output en- 
tropy and energy-constrained Holevo capacity and prove the 
additivity properties of two classes of Gaussian channels — 
the classical noise channels and thermal noise channels of ar- 
bitrary modes. 

The paper is organized as follows. In Sec.|lljwe introduce 
the notation and present basic facts about Gaussian states and 
the symplectic transformations used in this paper In Sec.lIIII 
we define Gaussian channels and introduce three figures of 
merits to quantify Gaussian channels — the maximal output 
p norm, the minimal output entropy, and the Holevo capac- 
ity. In Sec. II VI we formulate the additivity and multiplicativ- 
ity problems of Gaussian channels for Gaussian inputs. In 
Sec.|V]we prove a trace formula for symplectic eigenvalues 
and a majorization relation on symplectic eigenvalues that is 
an immediate consequence of the trace formula. By virtue 
of this majorization relation on symplectic eigenvalues, we 
prove the additivity and multiplicativity properties of classical 
noise channels and thermal noise channels of arbitrary modes 
in Sec. lVIl Section|yil|is devoted to concluding remarks. 



II. GAUSSIAN STATES 

In this section, we introduce the notation and summarize 
the basic facts about Gaussian states and symplectic transfor- 
mations 1 18]. We consider an n mode quantum system, such 
as a radiation field. Each mode corresponds to a quantum me- 
chanical harmonic oscillator with two canonical degrees of 
freedom and the quadratures of each mode correspond to the 
position and momentum of the harmonic oscillator. Thus an 
n mode state has a 2n canonical degrees of freedom. Let Qk 
and Pk denote the "position" and "momentum" operators as- 
sociated with the fcth mode {k - 1, 2, . . . , n). These operators 
or canonical variables are written in terms of the creation and 



2 



annihilation operators of the mode; 

Qk 

and 



(1) 



(2) 



where Uk denotes the energy of the feth mode {fi - 1). Since 
[aj,ak\ = [a\,a[] = and [aj,a[] = 6jk, we have [Qj,Qk\ = 
= and [Qj,Pk\ = Defining 



R - (RiR2,...,R2n) 



these canonical commutation relations (CCRs) can be written 
as [Rj,Rk] - i{Jn)jk- Here, J„ - ®"^^J\ with 



-l/2p ..T- 
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In the following, the characteristic function defined as xi^ - 
Tr[p'W(^)] plays a key role. Here, 'W(^) = exp(i^^J„R) is 
called the Weyl operators and p denotes the density opeartor 
The density operator in turn can be written in terms of its char- 
acteristic function and Weyl operators as follows. 



(5) 



A Gaussian state is defined as a state whose characteristic 
function is a Gaussian function: 



(6) 



Here, F > is a real symmetric matrix and D e R^". The first 
moment is also called the displacement or mean and given by 
m j - Tr{pRj) and the second moment is given by 

jjk = 2Tr\p{Rj -mj){Rk - ruk)] - i(J„)jk, (7) 

which is called the covariance of canonical variables. The 
2n X 2« real symmetric matrix {yjk) is called the covariance 
matrix y. T and D in Eq. (|6j are given by F = Jj^jJn and 
D — J„m. 

Note that due to our choice of canonical varibles, Rij-x - 

1/2 —1/2 

0)- Qj and R2j - a>j Pj (j - 1,2,..., n), the trace of the 
principal submatrix of the jth mode of y, 

gives the energy of the jth mode if in2j-\ = m2j - 0; 

^(^k (r2j-i,2;-i + r2;;2j) = Mk {{a]aj) + . (9) 

A density operator p is a positive semidefinite operator (p > 
0) with Trp = 1. The necessary and sufficient condition for 



p > in a Gaussian state is given in terms of the covariance 
matrix as follows uT 



y + iJn > 0. 



(10) 



Furthermore, the necessary and sufficient condition for a pure 
Gaussian state is given by f2^ 



dety = 1. 



(11) 



A linear transformation on canonical variables is written as 
R ^ R' - SR. Since the new variables R' also must conserve 
the CCR [R'j,R'i,] = iiJ„)jk, SJ„S'^ ^ J„ must hold. Such an 
2n X 2n real matrix satisfying S J„S ^ - J„ is called a sym- 
plectic transformation, S e Sp(2n, R) - {S\SJ„S^ - J„}, 
which forms a group so that 5 ' and 5 1^2 are symplectic if 
S,S 1,82 e Sp(2n, R). Furthermore, is also symplectic and 
det5 = 1 |21] if5 e Sp(2n, R). In this paper, we repeatedly 
use the following Williamson theorem 1221 . For a real sym- 
metric positive definite 2n x 2n matrix A = > 0, there 



(4) exists a symplectic transformation S e Sp(2n, R) such that 



SAS^ - diag(vi, vi, V2, V2, . . . , v„v„) 



(12) 



with Vj(> 0) called symplectic eigenvalues of A (j = 
1, 2, . . . , n). Equation il2\ is called the Williamson standard 
form of A. The symplectic eigenvalues can be computed via 
the eigenvalues of J„A, which are +ivj (j - 1, 2, . . . , n). Any 
symplectic transformation S e Sp(2n, R) can be decomposed 
into 

S = T^:^)ZT^^^ (13) 

with T^^\ G Sp(2n, R) n 0(2n) = K(n) and 

Z = diag(zi,zri,...,z„,z„i), (14) 

where Zj > 1 (j - 1,2, ...,n) |23]. 0(2n) denotes the or- 
thogonal group whose elements are 2n x 2n real orthogonal 
matrices. Equation jl3l l is called ffie Euler decomposition 
of symplectic transformations. K(n) is a maximal compact 
subgroup of Sp(2n,R) and is isomorphic to U(n), the unitary 
group whose elements are n x n unitary matrices |24]. The 
isomorphism is established via the following correspondance: 



T2i-\,lk-\ - T2j,2k - ReMj.i- - Uji^ 



and 



T2J- 



\,2k 



-T2j,2k- 



Imu 



'j,k^''j,k' 



(15) 



(16) 



where T € K{n) and uj^it are (j, k) components of n x n uni- 
tary matrices U. Using Eqs. ( I15I I and ( II6I 1. the isomorphism 
K{n) ^ U(n) is easily verified by direct calculations. 

Since the covariance matrix of an n-mode state is a 2« x 2n 
real symmetric positive-definite matrix, it can be cast into the 
Williamson standard form. In terms of symplectic eigenval- 
ues, condition jlOt is rephrased as Vj > 1 (j = 1, 2, . . . , «) and 
condition il l\ is written as vj - 1 for all j. 
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A canonical linear transformation corresponds to a unitary 
transformation in the Hilbert space. Such a unitary transfor- 
mation is defined by Us"W(^)Ul = "WiS "'^), and the density 
operator is transformed as p i-> UspU^^ - p correspondingly. 
It is easy to see that U^^ - U ^-i . The characteristic function of 
the new state pis given by Tr[t/spt/^W(^)] = Tr[p'W(5^)] = 
XiS^)- Accordingly, the covariance matrix and the displace- 
ment are transformed as y i-» S^^y(S^^)^ and m i-> S^hn. 
Note that the symplectic eigenvalues are invariant under such 
symplectic transformations on the covariance matrix. 

Coherent states, squeezed states, and thermal states are typ- 
ical Gaussian states, while the number states (of the single 
mode) given by \k} {k\ with 



\k}^-^(a^'f 10), 



k = 1,2,..., 



(17) 



are not. However, a vacuum state |0) {0|, which is a special 
case of the number states, is a Gaussian state with the covari- 
ance matrix. 



7v, 



diag(l,l). 



(18) 



This is the minimal-energy pure state. The coherent state is 
the displaced vacuum state so that the covariance matrix is 
given by Eq. jl8> but has a finite displacement. The thermal 
state of the single mode. 



i + <«)6^li + <«)j 



(19) 



(20) 



has the covariance matrix 

=diag(2 <«)-!- 1,2 <«) + l) 
with (n) being the mean photon number of the mode. 



III. GAUSSIAN CHANNELS AND THEIR 
QUANTIFICATION 

A Gaussian channel O is a completely positive trace pre- 
serving map that maps Gaussian input states p to Gaussian 
output states (S>(p) 1.25. .26.1 . The covariance matrix is trans- 
formed according to 



(21) 



where X and Y are 2n x 2n real matrices and Y is positive 
and symmetric (Y = Y^ > 0). The complete positivity of the 
channel is expressed in terms of these matrices as ll27ll 



Y + iJ„ - iX^J„X > 0. 



(22) 



Hereafter, we write a Gaussain channel by a capital greek 
letter and the coresponding transformation on the covariance 
matrix by the corresponding lower case greek letter. 

There are several figures of merits for quantifying quantum 
channels. Here we take three of them; the maximal output p 



norm '281, the minimal output entropy [29], and the Holevo 
capacity |30| for Gaussian state inputs. 

The Gaussian maximal output p norm is defined as 



^p(cD) = sup||(D(p)||^, 



(23) 



where ||p||^ = (Tr|p|'')"^ is the Schatten p norm (p > 1) with 

\A\ - VA' A. In Eq. (l23t . denotes the set of all Gaussian 
states. For a Gaussian state with covariance matrix y, 



Trp^ = n 



2P 



2'" 



FpivY 



where 



fpix) = (x+ir-(x-i)p. 



(24) 



(25) 



This formula has been originally derived in I'sT). Note that 
Trp^ is independent of the displacement m. We can verify that 
In/p is increasing and concave (Appendix A), so that Fp(v) - 
Yl"j=\ fpi'^j) is increasing and Schur-concave (Appendix B). In 
Eq. ( I23II. C? can be replaced by 0p, the set of all pure Gaussian 



states 



In terms of F,,, we have 



2" 



infFJvUiyp))]. 



The Gaussian minimal output entropy is defined as 
5^„(<D) = inf5((l>(p)), 



(26) 



(27) 



where S(p) - -Trplnp is the von Neumann entropy. Fol- 
lowed by the arguments presented in flTil . it can be shown 
that Q in Eq. j27t can be replaced by Qp. Since 



lim \\p\\^-S{p\ 
P^i+ dp ' 



(28) 



5tnin(<l>) can be computed through ^p(<l>). Note that 5'niin('l*) is 
also independent of the displacement m. Hereafter, we have 
occasions to write the von Neumann entropy as S (y) instead 
of S(py) when we are dealing with Gaussian states. 

By definition, the Holevo capacity for Gaussian state inputs 
or the Gaussian Holevo capacity is written as 



Cg(3>) = sup 



P-P(y.«,) L 



Smp))- j ^{dy,dm)S (0(P(y,,„))) 



where 



p(dy, dm)p( 



y,m) 



(29) 



(30) 



is the averaged signal state. In Eq. ( I29> . the supremum is taken 
over all possible probability measures p and signal states p(y,m) 
constituting the signal ensemble. Since the states p(y,„,) are in- 
finite dimensional states, the right-hand side of Eq. j29> be- 
comes any large number if we do not impose some constraint 
on the signal states. Here, we take the energy constraint. 



Trp// = & 



(31) 
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with 

// = J o., lala, + ^) = ^ 1^ (^2.-1 + 4,) ■ (32) 

k=i ^ z ^.^j 

Here we recall that the von Neumann entropy of a Gaussian 
state depends only on the covariance matrix, and that channel 
O affects only the covariance matrix. Therefore, if we find a 
single state P(y«,m) that minimizes S (<l*(p)), all possible Gaus- 
sian states p with the covariance matrix, y*, also minimizes 
S{^(p)). This observation indicates that the optimal signal 
ensemble that attains the Gaussian Holevo capacity consists 
of Gaussian states with the common covariance matrix y* and 
a certain probability distribution of the displacement m. If we 
restrict the signal ensemble to that described above, it suffices 
to take a Gaussian probability distribution for the probality 
measure fJ.{dy, dm) = fj.(dm). This is shown as follows I33il . If 
fi(dm) is a Gaussian distribution; 

u(dm) = e.xp(-m^ Y^^ m)dm. (33) 

with > 0, the averaged input signal state is calculated as 

P = J Hidm)p^^y^m) = P(y+Y,.0)- (34) 

That is, p is also a Gaussian state with the covariance matrix 
y - y + Y^ and has the vanishing displacement. Equation ( I34> 
even holds for > 0. Since the displacement of p is zero, 

TrpH ^\YjC^k ) + {rI)) ^^Ti-rw , (35) 

where y^^^j denotes the principal submatrix of the kl\\ mode 
of y defined by Eq. (|8|i so that the energy constraint [ QIH is 
written as 

n 

^w,Tr(r+i;)w =4£. (36) 

Since such a signal ensemble described above is not always 
optimal, we have 

Cg(0,£)> sup \sU{y+Y,))-S{cl,{y))\, (37) 

where the supremum is taken under the constraint ( I36> so that 
the right-hand side of Eq. is written as 

sup5(0[r]) - inf 5(0[r]) = sup5(0[r]) - 5,ni„(0), (38) 

y y y 

where the supremum is taken under the constraint 

n 

^ Trw,r[t] = 4£. (39) 

k=\ 

Here, we note the following extremal property of Gaussian 
states. For a given covariance matrix y, the von Neumann en- 
tropy is maximized for the Gaussian state L3 3^.34 1. Therefore, 



for the Gaussian signal state P(y,,„) and the probability measure 
p(c/y, dm), the quantity within the brackets of the right-hand 
side of Eq. \29\ cannot exceed the value of the right-hand side 
of J37I I. Therefore, the equality holds in the inequality J37> : 

Cg{<^,&) = sup 5 (0(r)) - S^,m. (40) 

y 

Again, the supremum is taken under the constraint [Eq. j39H. 

In the above arguments, we have assumed implicitly that 
there exists a Gaussian state with the covariance matrix satis- 
fying Eq. j39> . Otherwise, Cg(<1',£) should be zero because 
£ in the energy constraint j31> would be smaller than the sum 
of the energies of the zero-point oscillations over the modes. 



IV. ADDITIVITY AND MULTIPLICATIVITY PROBLEMS 
OF GAUSSIAN CHANNELS 

For the tensor product of the Gaussian channels, <1) — 
^"j^i ^j, it is evident from the definition that 

m 

^p(<i>)>n^''(*^^ (41) 

If the equality holds in the inequality ( I4H . we say that the 
maximal output p norm is multiplicative for Gaussian chan- 
nels Oj. To show the multiplicativity of the maximal output p 
norm, it suffices to show 

m 

inf Fp [v (0(rp))] = W inf Fp [v (^/T/,))] ■ (42) 

'' j=i '' 

In Eq. ( I42t . yp in the left-hand side of the equation is the 
covariance matrix of a pure Gaussian state on the composite 
Hilbert space "Hi ® ■ ■ ■ ®'H,„, while yp in the right-hand side of 
the equation is the covariance matrix of a pure Gaussian state 
on the Hilbert space 'Kj. 

By noting Eq. j28L it follows that if the maximal output 
p norm is multiplicative, then the minimal output entropy is 
additive. 

Let fi be the value for the energy constraint [Eq. j31H for 
the Gaussian channel and fi = '^"j=\ ^J- Fro^i the def- 
inition, the Gaussian Holevo capacity of the tensor product 
channel is greater than or equal to the supremum of the sum 
of the Gaussian Holevo capacity of individual channels; 

m 

Cg(<5,6)> sup yCc(<^j,6j). (43) 

Here, the supremum is taken over all possible combinations of 
Sj under the constraint Yj'j= \ - If the equality holds in 
the inequality J43l l. we say that the energy-constrained Gaus- 
sian Holevo capacity is additive for Gaussian channels O^. 
Now let p be a Gaussian state on the composite Hilbert space 
■Hi ® ■ ■ ■ ® Hn and define pj = 'Tr'^Hi<s: <s-Hj-imj^i» -e,n.fi-_^ 
noting the subadditivity of the von Neumann entropy 13511 . 
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S(P)< Z%iS(pj), we have 



S{c(,{y))<Y,s{^j(7j)), 



(44) 



where yj denotes the covariance matrix of the Gaussian state 
pj. Therefore, if the minimal output entropy 5'min(*l*) is addi- 
tive for the channels Oj, then 



Cg(0,£)< sup yCc(<J>j,&j). (45) 

(£.,KZ';L,fi,=£^ 

This implies the additivity of the energy-constrained Gaussian 
Holevo capacity 



Cg(0,£)= sup yCc(cD,,£,). 



(46) 



Serafini et al. proved that the Gaussian maximal output 
p norm of a tensor product of identical single mode Gaussian 
channels and that of single mode channels described by X, and 
y, \ M\V \ such that detX, are identical and F, > for all /, are 
multiplicative under Gaussian state inputs for p > \. Con- 
sequently, the Gaussian minimal output entropy and energy- 
constrained Gaussian Holevo capacity are additive for such 
tensor product channels. 



<.l-\,2m~\ 



in the Euler decomposition form [Eq. J13H to obtain 

2k k n 

j=l j=l l,m=\ 

k n 

+ Xi X ^2/-l,2m^'^m'22Z-l,2m 
7=1 /,m=l 
k n 

j=\ /,/«=! 
k n 

+ J]J]P'i,2n,^j'z-jQ2,,2m, (49) 



j=l /,/;,= I 

where zj > 1 (j - 1,2, ... , n). 



pU) _ 7^(1) 7^(1) , 7^(1)7^(1) 
^ ;,m -'2 j- 1 ,/-' 2j- 1 ,m 2 j/ 2 jjn 



(50) 



and 



with r(i\ e ^:(n). Using Eqs. O and (O for r*" and 
r^^\ the elements of f and Q are computed through the el- 
ements of « X n unitary matrices U and V as follows: 



2;- 1,2m 



(53) 



V. A MAJORIZATION RELATION ON SYMPLECTIC 
EIGENVALUES 



Lemma 1. - Let A be a 2« X 2« real symmetric positive- 
definite matrix (A = > 0). Then 



mm 

SJ„S 



Tr5 A5 ^ = 2 ^ vJ(A), I <k<n. (47) 

7=1 



The minimum in Eq. (|^) is taken over all 2kx2n real matrices 
S satisfying 5y„5^ = 7^. 

Proof. First of all, we note that a 2A:x2n matrix S satisfying 
S JnS ^ - Jkis the first 2k rows of a symplectic transformation, 
S',, e Sp(2n, R), and that A can be written in the Williamson 
standard form to obtain 



2k 



nun 

SJ„S'^=J, 



Tr5A5^= min V (5„£)a5I),j, (48) 

S„J„Sj,=J„ ^ 

7=1 



where we have used the fact that a product of symplectic trans- 
formations is a symplectic transformation and have defined 
Da = diag [v\(A), v\(A), vl(A), vl(A)]. Here, we write S„ 



Q2l-l,2.n-l = G2/,2,« = 2 "l^^^ i't^'n-P + "lA.p) ' ^^4) 



and 



e2/-1.2,„ = -e2;,2«,-l = "l^^^ hl'''>n,P + "'A) ■ (55) 



Substituting Eqs. (|53, (|54|i, and ^ into Eq. ^ yields 



2k 



7=1 

^ k n n 

" 4 X X ">''">,m(^' + Zr')(^m + Z;/) ^ v],(A)v,,^v;_^ 

7=1 /,m=l p=l 

+ 4 Zl «7;/«7;m(Z' -H Z^')(Z„ + Z„') ^ 

7=1 ',m=l ;'=1 

j=l Z,m=l ;j=l 



j=l /,m=l 
'=1 7=1 



P=l 



(56) 
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where * means complex conjugate. In the right-hand side of 
Eq. = t/Z+y, = WZ+V, W^^' = t/Z_y*, and 

^(4) ^ u*z_V, with 



Theorem 1. - Let A and 5 be 2nx2n real positive symmetric 
matrices (A = > 0, B = > 0). Then 



(57) 



Note that matrices W^'^D^W^'^^ are positive semidefinite; 

> (Z = 1, ... ,4). Now let be 

eigenvalues of Hermitian matrices . By the Schur 

theorem (Appendix B), we have 

k k 

2(W®Z)aH'«1)} . > ^Yw^'^DaW'^'^^ / = 1, . . . , 4. 

(58) 

Thus, we obtain 



v(A + B) v(A) + v(B). 
Proof. By Lemma 1 , we have 

k 

2Yv\A + B) = min Tr5(A + B)5^ 



(65) 



> min Tr5A5^+ min TrSBS^ 

SJ„S'^=Ji, SJ„S^=Jk 
k k 

= 2Y^v]{A) + 2Y^v]{B) 

= 2Y^[y{A) + ym]. (66) 



2*: 



4 it 



./=1 



4 /t 



By definition of the weak supermajorization r jB2l i1. inequality 
\66\ yields the desired relation \65\ . □ 



'=1 j=i 



Here, 



i](W<'>DAlV*'>"'") = a]{Z^VDaV^Z^) = ^J(Z+CZ+) (60) 

with C - VDaV^ > 0. The first equality is due to the unitary 
invariance of the eigenvalues of the Hermitian matrices. The 
eigenvalues, /lj(Z+CZ+), admit the following max-min repre- 
sentation isfill : 



VI. ADDITIVITY AND MULTIPLICATIVITY 
PROPERTIES OF GAUSSIAN CHANNELS 

In this section, we focus on two classes of Gaussian chan- 
nels; the classical noise channel and the thermal noise chan- 
nel. Both are important cases of Gaussian channels. 

In the classical noise channel, a classical Gaussian noise is 
added to the input states. Since 'W{^)Rj'WH^) = Rj + the 
classical noise channel is described by 



A {Z+CZ+) - max min 

^ wi,...,wj_ieC" A-#o 

A'_LVVi ,...,Vl',_ 



If we write y - Z+x, 



— —2 — • ^^^> n" Vdet Y J 



— P(y+Y.m) 



(67) 



\\x\\l^Yj(Z-')j\yj\' <J]\yj\' - Ml 



with F > 0. Namely, the transformations of covariance matrix 
is given by 



(62) 



>=i ./=i 



fir) -j + Y. 



(68) 



Hence, 



A.(Z+CZ+) > max min 

■' H'i,...,vvj_ieC" y#0 



y'Cy 



(63) 



Similarly, 



y±Z^'wi,...,Z_^'wj- 

= A](C)^y](A) 

so that /IJ(W<"DaW(""'") > vJ(A). 
iJ(W(2)DAW®t) > vJ(A). Since ^J(H'P>Da W<^>^) > 
and /lj(iy(4)z)^^(4)t) > 0, we find 

2k k 

Yj.SnDASl)jj>2Y^v]{A). 



In the thermal noise channel, the signal Gaussian states 
interact with an environment that is in thermal equilibrium. 
This channel is modeled by bemsplitters that couple the in- 
put Gaussian state and the thermal reservoir Let aj and bj be 
annihilation operators of the yth mode of the singnal state p 
and the thermal state p,/, that acts as a thermal reservoir The 
action of the bemsplitter is described by the transformations, 
aj i-> cos OjOj + sinOjbj and bj t-^ -sinOjaj + cos Ojbj. Ac- 
cordingly, the corresponding symplectic transformation takes 
the form. 



(64) 



cos 0jl2 sin 6*^12 
- sin 0jl2 cos 0jl2 



(69) 



7=1 



The equality holds for 5„ = Iin e Sp(2n, R). This completes 
the proof. □ 



Therefore, the output Gaussian state has the covariance ma- 
trix, 



"(y)^TrrH[S-'(y®yrh)(S-Yl 



(70) 
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where S = S ^ ® S 2 ® ■ ■ ■ ® S „, y,h = ®%i(2{n)j + l)l2 
denotes the covariance matrix of the thermal state p,/, with 
{n}j being the averaged photon number of the jth mode, and 
Tr,/, describes the trace over the thermal state. Using Eq. ( I69> , 
the right-hand side of Eq. ilQl is calculated as 



where v(Y) = 0" j v(Yj) with 



where 



and 



f,"'(y) = x'^yX + Y 



Y=^(2{n)j + m-rij)h 



(71) 



(72) 



(73) 



with rij = cos^ fl^ being the transmittivity of the bemsplitter. 
At zero temperature ((«)y - 0), the thermal noise channel is 
reduced to the lossy or attenuation channel i25l^3l . 



A. Classical noise channels 

For the /y-mode classical noise channel O^', the covariance 
matrix 7 is transformed according to <;*y'(y) - y + Yj, with 
Yj > 0. The tensor product of O^' is also a classical noise 
channel, and the covariance matrix of the output is given by 

111 

<^''(r) = 0<(r) = r + i', (74) 

where y is the covariance matrix of the input Gaussian state 
and Y - ©"Lj Yj. Since Yj is not always strictly positive 
definite, we add el2/, to Yj (s > 0); Yj(e) - Yj + slii. > 
i j - 1 , . . . , m) so that we can apply the Williamson theorem 
to Yj{e). Here, we write Y(s) - ® J^i i'j(e)- By Williamson 
theorem, there exists S j e Sp(2/j, R) such that 

S jY{E)jS ] = diag (yf{E), yfis), . . . , yj^^g), y^^^e)) 

= Dy,{s). (75) 

Here, we write S - 5^ so that SY{s)S^ - 

0™ J Dy^{s) = Dyie). By Theorem 1, we have 

v[y+ Yis)] ^v[SyS^ + Dyie)] v{y) + v(DKe)) . (76) 
Since Fp is increasing and Schur-concave, i76i yields 

Fp [y(y + Y(s))] > F„ (y(y) + v (Dy(s))) . (77) 
Here we can take the limit e — > to obtain 

inf Fp [v(f'(y))] > inf F,, {v{yp) + y{Y)) , (78) 



v(Yj) = lim(/f(s)Jj\e), 



(79) 



The infimum in the right-hand side of ilHi is achieved for 
v{yp) - (1,1, ...,1,1) and the equality holds if SyS^ takes 
the Williamson standard form. Namely, for the covariance 
matrix yp such that SypS^ = diag(l, 1, . . . , 1, 1), 

m Ij 

MFp[y{cp^>(yp))] = f] ^/^^^ + 



= f]infF4v(0f(rp))]. (80) 



That is, the maximal output p norm is multiplicative. Conse- 
quently, the Gaussian minimal output entropy and the Gaus- 
sian Holevo capacity are additive. Note that S defined above 
is the direct sum of local symplectic transformations and 
diag(l, 1, . . . , 1, 1) is the covariance matrix of the pure separa- 
ble state so that the optimal - 5"Miag(l, 1)(5^)"' 
is a separable pure state. This obvervation also indicates the 
multiplicativity of the maximal output p norm. 



B. Thermal noise channels 

For the L-mode thermal noise channel O'/', the covariance 
is transformed according to 



<'(r) = (r) + Yj, 



where 



(81) 



F, = y^(2<«)i + l)(l-/7,)l2 (82) 



and 



<pf\y)^X'yX+Y, 



(83) 



with X = V^2, Y = 0^^;(1 - m)h and < < 
1 . The tensor product of O^'' is a Gaussian channel and the 
covariance matrix of the output state is given by 

m 

<p"i7) = <P'j\7) = <P''°\7) + Y, (84) 

M 

where cf>'°\y) = 0" ; 0f (y) and Y = 0;'!^ Yj. Again, we 
add shn to y (e > 0); Y(e) - Y + sh,, to ensure Y(b) > 
(n = 2"Li hy Accordingly, we write (p"'{y,e) = (p^°\y) + Y{_e). 
By Theorem 1 , we have 

V [0"'(r, e)] V (0<°\r)) + V {Y{s)) . (85) 



8 



Since Fp is increasing and Schur-concave, ( I85> yields 

Fp{v [(jy'Hy, s)]] > Fp [v {cf \y)) + v {Y{s))\ , (86) 
Here we can take the limit e — ^ to obtain 

inf Fp [y {cf>'\yp))] > inf Fp [v (^("'(rp)) + y(Y)] . (87) 

Since the channel 0*"*(y) is completely positive, the Gaus- 
sian state with the co variance matrix ^^^\yp) is a physical 
state so that (0*"*(y/,)) > 1 (;' = 1,2, ...,«). For y - 

diag(l, 1, . . . , 1, 1), Vj f<^*°'(y)) = 1 and the equality holds in 
(gT). Hence, 

m 

inf Fp[v{cl>"\jp))] = nn^^ti+2(i-/7*)<«),] 

>=I k=l 
m 

= f[infF4v(^f(rp))]. (88) 

>=i 

That is, the maximal output p norm is multiplicative. Conse- 
quently, the Gaussian minimal output entropy and the energy- 
constrained Gaussian Holevo capacity are additive. 

VII. CONCLUDING REMARKS 

We proved the multiplicativity of maximal output p norm 
of classical noise channels and that of thermal noise channels 
of arbitrary modes for all p > I under the assumption that 
the input signal states were Gaussian states. As a direct con- 
sequence, we also proved the additivity of the minimal out- 
put entropy and the energy-constrained Holevo capacity for 
those Gaussian channels under Gaussian inputs. A majoriza- 
tion relation on symplectic eigenvalues was of importance in 
the proof. 

At present, very little is known about the inequalities related 
to symplectic eigenvalues of real positive-definite matrices. 
Efforts to unveil such unknown relations would assist in the 
analysis of entropic quantities of Gaussian states and would 
also shed light on the properties of Gaussian state entangle- 
ment l37ll and secure communication via Gaussian channels 

m. 
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APPENDIX A: CONCAVITY OF ln/p(x) 

It is readily seen that fp(x) is concave for 1 < p < 2 and is 
convex for p > 2. Therefore, In fp(x) is concave for I < p < 



2. In order to show the concavity of ln/;,(x) for /? > 2, we 
examine the second derivative of ln /p(x); 

where gp(x) = 4p(x^- 1)''^^ + fp(x)fp-2(x)- For p > 2, we find 
that gp(x) > so that In fp(x)/dx^ < 0. That is. In fp(x) is 
concave for p>2. Thus, \n fp{x) is concave for all p > 1. 



APPENDIX B: MAJORIZATION AND SCHUR CONVEXITY 

In this Appendix, we present definitions and basic facts on 
majorization and the Schur convexity (concavity) used in this 
paper fW] . 

For vectors, jc - (x\,X2, . . . , x„) wd y - (3'i,3'2, ■ ■ ■ ,y«) 
{xj,yj e R), we write x < y if xj < yj (j - 1,2, ...,n). 
Let x^ - (jcj, jCj, . . . , xi) denote the decreasing rearrange- 
ment of X, where xj > > ■ ■ ■ > xj:,. Similarly, let 
x^ = (xj, X2, . . . , xl) denote the increasing rearrangement of 
X, where x j < x^ < ■ ■ ■ < xJ . 

We say that x is majorized by y and write x < y if 

k k 

YuA-Yjyp k^l,2,-.-,n (Bl) 

with the equality for k - n. 

We say that x is weakly submajorized (weakly superma- 
jorized) by y and write x {<")y if 

^(^)iyf\ ^=l,2,...,«. (B2) 

It is easy to see that x <yif and only if x <^^■ y and x y. 

A real-valued function / defined on R" is said to be increas- 
ing if X < y fix) < f(y), while / is said to be decreasing if 
-/ is increasing. 

A real-valued function / defined on R" is said to be Schur- 
convex if x < y /(x) < f(y) while / is said to be Schur- 
concave if -/ is Schur-convex. 

A real-valued function / defined on R" satisfied x <w (<"' 
)y /(x) < {>)f{y) if and only if / is increasing and Schur- 
convex (concave). 

Let g be a continuous and nonnegative function on R. Then, 
f(^) - YVj=i g(^j) is Schur-convex (concave) if and only if 
logg is convex (concave). 

An application of majorization theory to matrix analysis is 
the following Schur theorem |40]. Let A e M„(C) be an Her- 
mitian matrix. Let diag(A) denote the vector whose elements 
are the diagonal entries of A and A(A) the vector whose coor- 
dinates are eigenvalues of A. Then, diag(A) < A(A). 
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